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We demonstrate that a multi-peak FMR spectrum, with lines corresponding to resonance in 
different ferromagnetic regions of a heterogeneous thin-film sample, can collapse to a single-peak 
spectrum if there exists a particular field configuration, or the configuration of the external magnetic 
field with respect to the film surface, in which dH las /dM ff — within the region magnetically 
dominating in the sample. 



I. INTRODUCTION 

The ferromagnetic resonance (FMR) spectrum of thin films is known to evolve with the angular con- 
figuration of the applied magnetic field with respect to the film surface. In this evolution the multi-peak 
FMR spectrum often becomes a single-peak one in a certain angular configuration, to regain its multi-peak 
character beyond it. A careful analysis of this effect, in which a multi-peak FMR spectrum 'collapses' to 
a single-peak one, leads to the observation of two types of behavior of the collapsing spectrum. In one 
type the intensity of all the peaks except the first one diminishes progressively to vanish completely in 
the collapse configuration. Beyond this configuration angle the 'satellite' resonance lines emerge again in 
an unchanged order. The other type of collapse involves resonance positions rather than intensities: in 
the collapse configuration the lines shift to the position of the main peak, to reemerge in a reversed order 
beyond this configuration. 

The first type of collapse has been long known in the literature to be a surface effect due to the changes 
of the surface magnetic anisotropy (responsible for the surface spin pinning) with the configuration of 
the external field with respect to the surface of the thin film. In this pattern a complete collapse of the 
resonance spectrum occurs in a specific angular configuration (referred to as critical angle of the surface 
anisotropy) in which the surface anisotropy has no effect on the surface spin pinning, and the resonance 
precession of the all spins is homogeneous throughout the sample. In contrast, the other type of collapse is 
a bulk effect that occurs as a result of the ferromagnetic resonance in separate regions of slightly different 
magnetic character in a heterogeneous thin-film sample. A heterogeneous magnetic structure is known to 
result in heterogeneous conditions of ferromagnetic resonance in the sample. Thus, in this interpretation 
each line in the multi-peak resonance spectrum corresponds to a resonance in a different region of the 
sample, and the collapse configuration of the external field is a particular configuration in which all the 
magnetically different regions participate in the resonance for resonance fields that differ only negligibly 

The collapse of the FMR spectrum in this specific angular configuration of the external field in which 
the resonance heterogeneity of the sample is eliminated has not been thoroughly analyzed in the literature 
so far. The aim of this paper is to elucidate in detail the theoretical grounds of this effect. 

The paper is organized as follows. In Section II we recall the derivation of the universal Smit-Beljers- 
Suhl formula expressing the condition for ferromagnetic resonance to occur in a homogeneous bulk fer- 
romagnetic sample. On the basis of this formula, in Section III we derive the configuration condition 
of resonance in thin films. In Section IV we use this condition for analyzing the configuration evolution 
of the resonance field in thin-film samples with a perpendicular uniaxial anisotropy. We formulate a 
universal condition for a multi-peak FMR spectrum to collapse into a single-peak one in such samples, 
and (in Section V) demonstrate the bulk character of this effect. 



II. SIMPLIFIED DERIVATION OF THE SMIT-BELJERS-SUHL RESONANCE FORMULA 

Let us consider the dynamics of a body with an angular momentum J and a magnetic moment m = jL 
collinear to it (7 is the gyromagnetic ratio) in a magnetic field H . The field acts on the dipole to produce 
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a torque m x H that sets the body in motion according to the equation: 

rh = jm x H . 



(1) 



Let us represent the vectors in the above equation in the local basis of orthonormal spherical vectors 
r,9,4> determined by the vector m(m,9,<f)) (obviously, m = mr, see Fig.l and comments in [1]). Thus: 



rh = 7(mr) x (H r r + H g 6 + H^). 



(2) 



Now, let us determine the components Hg and of the magnetic field H by considering the change in 
the energy of the dipole with an infinitesimal rotation (we only consider rotation by the polar angle 8). 
The energy of a magnetic moment m in a magnetic field H is expressed by the equation: 



E(m) = rn ■ H, 
which, represented in the basis r, 0, </>, becomes: 



(3) 



E(m) = -(mr) • (H r r + H g 6 + Hrf). 
As a result of the rotation of m by a small angle Ad (m — > m! = mr') the energy changes to: 

E(m') = -(mr') ■ (H r f + H e 9 + H+jf). 
The vector r' can be expressed as f + A66, so the above equation will become: 

E{m') = -m(r + AO0) ■ (H r f + H e + H^). 
The following relations are seen to occur: 



AE = -mH e A6 



1 dE 
m~dO 



(4) 
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Now, we can get back to the equation of motion (1) and write the cross product as the determinant: 

(9) 
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Hence we obtain the system of equations: 



m 

-6 sin9 
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(10) 



0,<t> 



This system of equations is to be solved in four steps. We (i) determine the equilibrium angles 9q 
and 4>o as the solution of the system of equations dE/d9 — dE/dcf) = 0; (ii) expand E and sin6> into 
a Taylor series at 9q and (j>o, respectively; (iii) assume that the time dependence of the angles has the 
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form e lLOt , and (iv) exclude all the terms except the lowest-order ones (see [2]). This procedure leads to 
the Smit-Beljers-Suhl formula [3, 4] that represents the ferromagnetic resonance condition: 



u _ 1 
7 m sin 9 

where the values of the derivatives correspond to the equilibrium angles 9o and 0o- (Note by the way 
that the radicand is the Hessian of the function E at the equilibrium point; its positive value is indicative 
of the existence of a minimum of the function at the point at which the first partial derivatives vanish.) 
Since we shall henceforth consider samples with a magnetization M ', the formulas derived above, with M 
in place of m, will apply to the description of the magnetization dynamics in the studied case. 
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III. RESONANCE CONDITION IN THIN FILMS WITH UNIAXIAL ANISOTROPY 

Let us consider a sample in which the free energy density E is the sum of the Zeeman energy the 
demagnetization energy and the uniaxial anisotropy energy: 

E = M ■ H + 2tt(M -nf - K(M -u/M) 2 , (12) 

where n is a unit vector normal to the surface of the sample, and it is a unit vector oriented along the easy 
magnetization axis. The applied magnetic field and the magnetization of the sample will be henceforth 
represented as (see Fig. 2a): 

M = M(sin 9 cos 4>, sin 9 sin (/>, cos 9) 

H = H(0,sm6 H ,cos6 H ). ^ ' 

We shall consider three special cases corresponding to three different choices of the easy axis ([100], [010] 
or [001], see Fig. 2b). In each case we shall investigate the angular configuration dependence H ics (9h) 
implied by the assumed form of the free energy density (12). 



A. Perpendicular uniaxial anisotropy, u — [001] 

In this case the unit vector n normal to the surface of the sample and the unit vector u oriented along 
the easy axis are identical: 

n = u=[mi). (14) 

Thus, the free energy density reads: 

E = -HM (sin <f> sin 6 sin 6 H + cos 6 cos H ) + (2irM 2 - K) cos 2 9. (15) 

Let us determine the equilibrium conditions of the system. The first partial derivatives with respect to 
the angles 9 and 4> have the form: 

dE 

— = ~MH (sin <f> cos 9 sin 9 H - sin 9 cos 9 H ) - (2irM 2 - K) sin 29 

dE ^ 
— — = — MH cos sin # sin . 

d(j> 

From the condition of vanishing of these derivatives we obtain the equations that determine the coordi- 
nates 9m and <pM of the equilibrium point (Fig. 2b): 

2Hsm{9 M - H ) = 4irM cS sin26» M (17a) 
cos (pM = <= <j> M = tt/2, (17b) 



4 



where M c s is the effective magnetization, defined as: 

47rM cff = 4ttM - 2K/M. 
In this case the Smit-Bcljcrs-Suhl formula becomes: 



= [Hcos(9 M - 9 H ) - 47rM off cos 2 M ] [Hcos(6 M - H ) - 4nM eS cos 20 M ] ■ 



(18) 



(19) 



In order to plot the corresponding configuration dependence H tcs (0h) of the resonance field we must 
solve the system of two nonlinear equations (17a) and (19). 



B. In-plane uniaxial anisotropy, u — [010] 

In this case the unit vector n normal to the surface and the unit vector u oriented along the anisotropy 
direction are: 

n = [001] and u = [010]. (20) 

The free energy density is: 

E = -HM (sin (f> sin sin H + cos 9 cos 9 H ) + 2irM 2 cos 2 6-K sin 2 9 sin 2 (f>. (21) 
The first partial derivatives with respect to the angles 9 and <\> read: 



dE 
~89 

dE 



= — MH (sin <p cos 9 sin Oh — sin 9 cos Oh) 

- AttM 2 sin cos - 2K sin 2 sin 9 cos 9 
= — MH cos (f> sin 9 sin Oh — 2K sin 4> cos </> sin 2 0. 



Thus, the equilibrium direction of magnetization is determined by the equations: 

2K 

i0 H ) + (4ttM + — sin 2 (t> M ) sin 20 M 



= 2i7(sin <pM cos 0m sin 0// — sin#M cos ( 



( 2K 
sin M \H sin 6> ff + sin </> M sin # M 



The 



= COS 4>m 

Smit-Beljers-Suhl formula for this case is: 
Hcos(0m - Oh) 



- (A-KM + j COS 20 M 



' „ sin H 2K 
H . " + — 
sin &m 



or, equivalcntly: 

• 2 r / 2K' 

Hcos(9m-9 h )- U7rM+ — 



T 2 



cos 20m 
Hcoa(0 M 



x 



/ 2K\ 2K 

- H ) - [A-kM + — J cos 2 M + — 



(22) 



(23) 



(24) 



• (25) 



Note that this condition does not involve the azimuth angle <p, since the equilibrium condition (23) 
requires again that 4>m = 7r/2. 
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C. In-plane uniaxial anisotropy, u = [100] 



In this case the unit vector n normal to the sample surface and the unit vector u parallel to the 
anisotropy direction are: 

n = [001] and u = [100]. (26) 

The free energy density reads: 

E = —HM (sin 6 sin 9 sin 9 H + cos 9 cos 9 H ) + 2irM 2 cos 2 6 - K sin 2 9 cos 2 6. (27) 
The first partial derivatives of the energy density with respect to the angles 9 and 6 have the form: 

dE 
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— HM (sin 6 cos 9 sin 9h — sin 9 cos 9h) 

- AirM 2 sin 9 cos 9 - 2K cos 2 6 sin 6 cos ( 



(28) 



dE 2 
— — = — HM cos 6 sin 9 sin 9h + 2K sin 6 cos 6 sin 9. 
am 



Thus, the equilibrium direction of magnetization is determined by the angles 9m and S M that fulfill the 
equations: 



' = — HM (sin 4>m cos 9m sin 0// — sin 9m cos ) 
< — 47rM 2 sin 9m cos 9m — IK cos 2 4>m si n $m cos #m 
k = — HM cos ^>m sin 9m sin 0# + 2_ftT sin 4>m cos <^m sin 2 9m- 

The Smit-Bcljers-Suhl formula for this case has the form: 



(29) 



2K 

H (sin <j> M sin M sin 6* ff + cos 6m cos ) — ( A-kM + -— cos 2 0m ) cos 29m 
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cos 2(f) 



M 



H cos 6 m cos 6m sin 6h 2K 



or, equivalcntly: 



sin# 



M 



sin 24>m cos 6*m I (30) 



2K \ 

H (sin 0m sin sin 6*^ + cos 6m cos 6*^) — ( 4irM + cos 2 <f>M J cos 2^m 

/ 2K \ 2K 

H(sm 6m sin 9m sin + cos 6m cos ^ ) — ( 4ttM H cos 2 6 M ) cos 2 9m H cos 20m 

\ M J M 



2K 
~M 



sin 2 4> M cos 2 0m cos 2 6>m ■ (31) 



IV. THE EXISTENCE OF A RESONANCE INTERSECTION POINT IN 
HETEROGENEOUS THIN FILMS 



In this Section we shall only consider the simplest of the three cases mentioned above, i.e. a thin film 
with perpendicular uniaxial anisotropy (case [001]); we leave the remaining two cases to be similarly 
discussed in a separate paper. The condition (17a) allows to determine the equilibrium angle 9m of 
the magnetization vector corresponding to a specific configuration of the external field 9h- Figure 3 
shows the consequent configuration dependence 6m(6h) f° r two values of 47rM cf f. The deviation of the 
magnetization angle 6m from the field angle 6h is seen to be the largest in the middle of the interval 
[0, 7r/2] , and to grow with increasing value of 47rM e ff; in this range the deviation is of ca. 10° to 15°, far 
from being negligible! 
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Now, let us assume that the thin film under consideration is heterogeneous and planarly stratified into 
regions with two different values of effective magnetization AirM e g. The resonance condition (19) - with 
the equilibrium condition (17a) used - leads to two respective configuration resonance curves iJ res (^ff), 
see Fig. 4a. This means that for a configuration Oh each of the ferromagnetic strata resonates for a 
different value of the external field. This applies to the whole range of Oh with the exception of one 
particular configuration corresponding to the point of intersection of the two curves. In this particular 
configuration both regions - in spite of their significantly different magnetic parameters! - produce 
a resonance simultaneously for the same magnitude of the external field. This implies homogeneous 
precession throughout the sample, as if its heterogeneity were completely eliminated in this configuration! 
Thus, for this particular angle Oh the sample can be regarded as dynamically homogeneous. The dynamic 
homogeneity (DH) is evidenced by the collapse of the two-peak resonance spectrum to a single resonance 
line. Note that past this particular angle the two peaks reemerge, but the resonance sequence is reversed. 

The occurrence of a dynamic homogeneity angle (DHA) is an intrinsic characteristic of the Smit-Beljers- 
Suhl (SBS) resonance formula. Thus, we have grounds to ask whether the SBS equation can provide the 
basis for an analytical condition that would allow to predict the DHA angle for a ferromagnetic thin- 
film sample. We think we have managed to formulate such a condition a posteriori on the basis of the 
numerical results shown in Fig. 4b, presenting the configuration dependence of the derivative dH/dM e s 
determined for each of the ferromagnetic regions shown in Fig. 4a. The DHA angle is seen to be in 
the range delimited by the extreme configurations for which the considered derivative vanishes in each 
stratum. The range shrinks with decreasing difference between the extreme values of A-kM c r in the two 
regions. Thus, we can expect that the DHA angle will be determined with a good approximation by the 
condition of vanishing of the derivative: 

(32) 

in the region the magnetic properties of which predominate in the sample. 



V. CONCLUSION 



It is worthy of notice that the above-discussed collapse of the FMR spectrum in a particular field 
configuration defined by the DHA angle is a dynamic volume effect, since only bulk quantities figure in 
the adopted expression for free energy density, and the presence of the surface is only manifested in the 
shape anisotropy of the sample. Thus, the energy of surface anisotropy was not taken into account in 
our considerations. This is an important statement, since also the surface anisotropy by itself can cause 
a multipeak SWR spectrum to collapse to a single resonance line in a certain configuration, referred to 
as critical angle [5], of the applied field with respect to the film surface. However, there is a substantial 
difference between the collapse of an FMR spectrum and that of an SWR spectrum. In the former the 
resonance line positions merge to form a single line at the DHA angle, while the collapse of an SWR 
spectrum consists in the suppression of the intensity of all the lines except the main one at the critical 
angle. This is due to the different nature of the two effects: the essence of the FMR collapse is the 
attainment of homogeneous resonance dynamics throughout the volume of the sample, whereas the SWR 
effect consists in the elimination, in the critical field configuration, of the impact of the surface anisotropy 
on the spin dynamics in the thin film. Therefore, we shall refer to the latter configuration as surface 
critical angle (SCA). The SCA is sensitive to the conditions on the surface of the sample and can vary 
with them, while the DHA, only determined by the bulk characteristics that enter the condition (32), is 
a fixed parameter of the sample. It should be emphasized that the above-mentioned differences between 
the two effects can be used in practice as a criterion that allows to determine whether a collapse of a 
multipeak FMR spectrum observed in a thin film is of surface or bulk character. 
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Fig.1 



FIG. 1. Local basis of spherical vectors r, 0, (f> related to vector m. 




FIG. 2. Definition of the angular configuration of the magnetization M and the applied magnetic field H with 
respect to the film surface. 




EXTERNAL FIELD ANGLE, 6 U 



FIG. 3. Equilibrium magnetization angle 6m vs. external field angle Oh determined for a thin film with per- 
pendicular uniaxial anisotropy (f=33.5 GHz); two curves corresponding to two effective magnetization values are 
plotted for comparison. 
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FIG. 4. Configuration dependence of (a) the resonance field H ics (9h) and (b) its derivative dH las /dM a ff in a thin 
film with perpendicular uniaxial anisotropy. (a) Note the intersection of the resonance curves plotted for the two 
regions of different effective magnetization magnitude. The intersection point corresponds to the particular field 
configuration in which the sample produces a homogeneous resonance for the same field magnitude throughout 
its volume, (b) Note that the homogeneous resonance configuration lies in the range determined by the vanishing 
of the derivative dH res /dM e g in each resonating region. 
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